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Abstract 

Single ionization of N2, O2 and F2 by intense laser fields is simulated using a time-dependent density 
functional approach within the exchange-only local density approximation. In particular, the response to 
laser pulses having a wavelength of A = 390nm (corresponding to a frequency-doubled Ti:Sapphire laser) 
is investigated. Single ionization yields of N2 and F2 are seen to agree to within a factor of 3, in keeping 
with experimental results but in disagreement with predictions based on the symmetry of the highest occu- 
pied molecular orbital. The underlying mechanism for this appears to be the importance of multi-electron 
molecular effects in the single ionization of F2. In addition single ionization of O2 is calculated assuming 
the ground state configuration to be either a singlet or triplet state. It is discovered that the single ionization 
process is largely insensitive to the multiplicity of the initial state. 
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I. INTRODUCTION 



In the last two decades the interaction of atoms and molecules with intense laser pulses has 
attracted widespread research. This highly non-perturbative interaction results in a number of 
non-linear processes such as ionization, harmonic generation and above-threshold ionization. At 
the heart of all these processes is single electron ionization which gives rise, either directly or 
indirectly, to the other processes. 

Compared to atoms, molecules represent a more complex class of system due to their multi- 
centre nature which introduces additional vibrational and rotational degrees of freedom. It is per- 
haps surprising therefore that the ionization of molecules by intense laser pulses share many of the 
characteristics with the ionization of atoms. The mechanisms of ionization can be characterized 
as either multiphoton transitions or tunnelling ionization, or some combination of both. The pro- 
cess can be classified into either the tunnelling or multiphoton regime by the Keldysh parameter, 



defined as = v /p/2f/p, where the internal binding energy is Ip and the external laser-driven 
kinetic energy is Up. In the long- wavelength, intense field limit 7^ ^ 1 and tunnelling dominates 
and thus the excited states spectrum should play no role in the ionization process. 

In a range of early experiments this was found to be the case: it was discovered that single ion- 
ization rates for molecules are roughly identical to those of noble gas atoms provided the ionization 
potentials were comparable llll^lsl ^. In one of these experiments lO] the single ionization rates 
of O2 (Ip = 12.07eV) and of its companion noble gas atom xenon (Ip = 12.13eV) were measured 
simultaneously by exposing a target containing a mixture of xenon atoms and O2 to a laser pulse 
of wavelength A = 10.6/im. The ionization rates were found to be similar. 

However, in later experiments carried out at Ti:Sapphire laser wavelengths (A ~ 800nm) it was 
found that while most atom-molecules pairs obeyed this finding - such as N2 (Ip = 15.58eV) and its 
companion atom argon (Ip = 15.76eV) - a number did not. For example, the single ionization rate 
of O2 was suppressed with respect to ionization rate of xenon by several orders of magnitude jslQ] 
while ionization of D2 was also suppressed with respect to its companion noble gas atom, argon. 

Considerable interest has been generated by these findings and a number of explanations have 
been put forward for the origin of the suppression, most notably in O2. Talebpour et al o] sug- 
gested a dissociative recombination process leading to a decrease in the single ionization signal. 
However, later experiments lO] concluded that dissociative recombination cannot be the cause of 
suppression since the suppression was observed for a range of laser polarization ellipticities. An 



2 



n 

electronic correction to tunnelling theory was suggested by Guo [7] and reproduced the correct 
suppression of O2, provided that correct parameters for an effective nuclear charge and ionization 
potential are chosen. _ 

Muth-Bohm et al [8] explained the suppression of the O2 single ionization rate using a general- 
ization of intense-field many-body ^-matrix theory (IMST) which included an interference term. 
They showed with their calculations that the suppression of the O2 signal with respect to xenon, 
and its absence in N2/Ar, is due to a symmetry induced dynamical effect whereby interference be- 
tween ionizing wavepackets emitted from the two distinct nuclear centres is either destructive or 
constructive in the low energy limit, depending upon the symmetry of the highest occupied molec- 
ular orbital (HOMO). Thus for O2, in which the HOMO has an anti-bonding character (ground 
state configuration lcr|lcr^2cr^2(7^3(T^l7r^l7r|), single ionization should show suppression. For N2, 
in which the HOMO has a bonding character (ground state configuration lcr|lcr^2cr^2(T^l7r^3(Tg), 
single ionization should not be suppressed. In addition, according to the symmetry argument of 
IMST, when suppression does occur we can write 

kN-R<^l, (1) 

where k]^/2 = Ntu — Up — Ip is the kinetic energy of the electron on absorbing photons, oj the 
laser frequency and R the bond length of the molecule. Thus, if only the wavelength is varied, the 
suppression should be enhanced at shorter wavelengths. This explanation is consistent with the 
02/Xe results where suppression of O2 was observed at A = 800nm but not at A = 10.6/im. 

Following the symmetry argument it was postulated that since F2 (Ip = 15.69eV) has an ion- 
ization potential similar to N2 and argon but has valence electrons with the same symmetry as 
O2 (ground state configuration of F2: lcr|lcr^2(j|2cr^3(T|l7r^l7r^) ionization of F2 should also be 
suppressed with respect to either N2 or argon Isl]. However a later experimental study Isl showed 
that the ionization of F2 is not suppressed with respect to that of N2. 

In this paper we study the response in time of N2, O2 and F2 to laser pulses having a wavelength 
of 390nm and 300nm. We find single ionization suppression in O2 and its absence in F2, in ac- 
cordance with the experimental results at A = 800nm. Within our framework of a time-dependent 
density functional approach we are able to explain the deviations from the predictions based on 
the symmetry of the HOMO. The paper is arranged as follows. In section |ll| the time-dependent 
density functional approach is set out and the procedure for calculating single ionization rates de- 
scribed. In section |in| the time-evolution will be considered. Finally, some conclusions will be 
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drawn in section HVl 



n. METHOD 

The time-dependent density functional method provides the most detailed, practical and feasi- 
ble ab initio approach for tackling many-body problems. Density functional theory (DFT), as first 
introduced by Hohenberg and Kohn llflll and Kohn and Sham [11] describes a system of interacting 
particles in terms of its density. The theory is based on the existence of an exact mapping between 
densities and external potentials and leads to the density of the interacting system being obtained 
from the density of an auxiliary system of non-interacting particles moving in an effective local 
single particle potential, i.e. the particle interactions are treated in an averaged-over manner. A 
time-dependent formalism of DFT (TDDFT) was provided by Runge and Gross 13], who showed 
that the time-dependent density could be obtained from the response of non-interacting particles to 
the time-dependent local effective potential. In principle, many-body effects are included exactly 
through an exchange-correlation functional; in practice, the form of this functional is unknown 
and at best it can only be approximated. 

Such an approach has been widely used in treating the interaction of molecules with intense, 
short-duration laser pulses. For instance, in the approach of Chu and co-workers lisifl^ a pseudo- 
spectral mesh technique is used in the solution of the Kohn-Sham equations. A non-adiabatic 
quantum molecular dynamics (NA-QMD) method, in which the electron dynamics are treated 
quantum mechanically using TDDFT, has been developed by Schmidt and co-workers jljs]. In this 
work the electron orbitals are expanded in terms of Gaussian basis functions. More recently, Castro 
et al [ 17] developed grid-based NA-QMD technique to study harmonic generation in H2. However, 
these calculations were limited to a ID treatment of the electron dynamics. Another grid-based 
DFT method has been developed by Otobe et al lisl . While this method has calculated tunnel 
ionization rates for N2, O2 and F2 two approximations have been made. Firstly, a pseudo-potential 
description of the electron-nuclear potential was used meaning that only the valence electrons in 
the molecules were treated. Secondly, only static field ionization rates were calculated. 

The current approach utilizes a NA-QMD approach in which the electron dynamics are de- 
scribed by a hybrid finite difference-Lagrange mesh technique lisl . For the present calculations, 
however, a fixed nuclei description has been employed. Details of the approach are now given. 
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A. Time-dependent Kohn-Sham equations 



In the TDDFT approach, the total Ai'e-electron Kohn-Sham wavefunction is written as a single 
determinant of one-particle Kohn-Sham orbitals. Denoting the spin state of each orbital by the 
label a =t, |, we can write the electron density as 

Nu 

n{r,t)= Y,Mr,t)= |^..(r, (2) 

where is the number of electrons in spin state a and 'ipiai'f, t) are the Kohn-Sham orbitals. The 
orbitals are obtained through the solution of the time-dependent Kohn-Sham equations 



•5 , . ^ 



-^V2 + Kff(r,t) 



i^ia{r,t), (3) 



where 

V,sir, t) = Kxt(r', R, t) + Vnir, t) + V^c.{r, t), (4) 
is the time-dependent effective potential which is given in terms of the external potential 

Ven{ri, R, t) = Vion.iVi, R, t) + f/eleclT"*, t) , (5) 



where f/eiec {ri-, t) denotes the interaction between electron i and the applied laser field and where 

\Ri - r 



l^onslr-i, R,t) = Y, VioniVi, Ri, t) = -J^ | „ y 



1=1 1=1 

denotes the Coulomb interaction between electron i and all ions. The time-dependent effective 
potential also depends on the Hartree potential 

VH{r,t)= [dr'p^^, (7) 
J \i — ''"I 

and the exchange-correlation potential 

(^K_„ fr),+ . nil 

(8) 



6n^ 

where E^c [n^, raj is the exchange-correlation action. 



na=na{r,t) 



B. Treatment of the exchange-correlation potential 

All many-body effects are included within the exchange-correlation potential, which in practice 
must be agproximated. While many sophisticated approximations to this potential have been de- 
veloped pi9l], the simplest is the adiabatic local density approximation in the exchange-only limit 



(xLDA). In this case the exchange energy functional is given by 

1/3 



(9) 



from which the exchange-correlation potential 

gXl/3 



V.cAr,t) = -{-) nl/%r,t), (10) 



TT 

can be obtained. This approximate functional is easy to implement and is one of the most widely 
used exchange-correlation functionals. However, it does suffer from a number of drawbacks, most 
notably it contains self-interaction errors. This self-interaction means that the asymptotic form 
of the potential is exponential instead of Coulombic. Therefore, the electronic properties and 
response of the system can differ markedly from those of the actual system (as we shall see, for 



example in section UlI All . 



C. Numerical details 

Precise numerical details of how the code is implemented are given in IirI . Briefly, the nu- 
merical implementation of TDDFT uses a cylindrical grid treatment of the electronic Kohn-Sham 
orbitals. As in 11221] a finite difference treatment of the ^-coordinate and a Lagrange mesh treatment 
of the p-coordinate based upon Laguerre polynomials is used. For the case of diatomic molecules 
and considering a linearly polarized laser pulse with the laser polarization direction parallel to 
the molecular axis the azimuthal angle can be treated analytically. The time-dependent Kohn- 
Sham equations of TDDFT are discretized in space using these grid techniques and the resulting 
computer code parallelized to run on massively-parallel processors. Several parameters in the 
code affect the accuracy of the method. These are the number of points in the finite difference 
grid (A^^), the finite difference grid spacing (Az), the number of Lagrange-Laguerre mesh points 
(Np), the scaling parameter of the Lagrange-Laguerre mesh (hp), the order of the time propaga- 
tor (Nt) and the time spacing (At) llisll . Li all the calculations presented here, converged results 
were obtained using the following parameters (atomic units are used throughout): A^^ = 2291, 
Az = 0.05, Np = 43, hp = 0.28838771, Nt = 18 and At = 0.02. The code was parallelized and 
the calculations were carried out using 79 processors. 
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D. Determination of Ionization 



Within TDDFT all observables are functionals of the electronic density. However, as in the case 
of the exchange correlation functional, the exact forms of these functionals may not be known and 
must therefore be approximated. The functional describing ionization falls into this category. To 
date, most calculations of ionization within TDDFT have been obtained using geometric properties 
of the time-dependent Kohn-Sham orbitals Briefly, an analysing box is introduced as a way 
to approximately separate the bound- and continuum-state parts of the wavefunction. We define 
the box such that all relevant bound states of the wavefunction are contained inside this box while 
the continuum contributions are found outside the box. In that case the number of bound electrons 
is given by 

A^bound(t) = / dr n{r,t), (11) 

J inside box 

while the number of continuum electrons is given by 

A^escW = / drn{r,t). (12) 

J outside box 

We then define the number of bound and continuum electrons for each Kohn-Sham orbital as 
(dropping spin subscripts) 

Njit)= [ dr mr,t)\\ (13) 



inside box 



for bound electrons and 



iV,(t)= / dr mr,t)\\ (14) 



outside box 



for continuum electrons. Thus we obtain approximate ion probabilities, P (t), for charge state k. 
In particular we find 

P%t)=N,{t)---Nr,M, (15) 



and 



P\t)=Y,Ni{t)---Nn{t)---NNM- (16) 



n=l 



III. RESULTS 

We now study the response of N2, O2 and F2 to intense laser pulses. Firstly, we compute 
the ionization potentials of the various molecules. Then we compare single ionization of N2 and 
F2 showing that the F2 yield shows no suppression with respect to the N2 yield. Thirdly, we 
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compare the yields of O2 using two initial configurations of the molecule, namely a singlet and 
a triplet state. This shows that the suppression of the signal is not due to the fact that the O2 
ground state is a triplet configuration. In order to gain a better understanding of the electronic 
dynamics we analyse the time-evolution of the Kohn-Sham orbitals, showing that in the case of F2 
the response of the core electrons show a significant response to the pulse. 



A. Ionization Potentials 



The starting point for any simulation of the response of diatomic molecules to intense laser 
pulses is an appropriate description of the field-free structure. In our calculations we assume 
that the molecules are initially in their ground states. Starting from the equilibrium ion separa- 
tion the ground state electronic density is calculated self-consistently from the time-dependent 



Kohn-Sham equations using an iterative Lanzcos method ll2lL 22! . Of particular importance in the 



present calculations are the ionization potentials of the different molecules. These are obtained by 
first calculating the ground state energy of the particular molecule and then that of the molecule 
with one electron removed. The results are presented in tableUfor N2, O2 and F2. We see that the 
ionization potentials of N2 and O2 compare well to the experimental values ll23h while the ioniza- 
tion potential of F2 shows a 10% difference with the experimental value. The difference can be 
explained in terms of the choice of exchange-correlation potential, as discussed earlier. Obviously 
in comparing single ionization in N2 and F2 the errors in their ionization potentials will have an 
important impact upon their subsequent response. It has been pointed out to us I24I1 that the 10% 
decrease in the ionization potential of F2 could increase the single ionization yield by two orders 
of magnitude. However, such an estimation, obtained using IMST, only takes into account the 
response of the HOMO to the field. As our results will show, it appears that molecular effects due 
to the other orbitals are extremely important in this system and such a large increase in ionization 
yield due to a lowering of the ionization potential is unlikely. This is supported by the static field 
tunnel ionization calculations of Otobe et al |15] who employed a DFT approach using the self- 
interaction free KLI approximation. In these calculation the ionization potential of F2 was much 
closer to the experimental value. However, no suppression of the F^ signal was evident. 
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B. Single ionization yields of N2 and F2 



TableHUpresents a comparison of the single ionization yields for and after the interaction 
of the neutral molecules with a 24 cycle laser pulse (pulse length, r = 3 1 .2fs) having a wavelength 
of A = 390nm over the intensity range / = 1 x lO^'* W/cm^ to / = 8 x 10^'* W/cm^. As discussed 
in the introduction, if the symmetry induced dynamical effect is responsible for ionization sup- 
pression, equation (1) predicts that the yield of F2 ions should be suppressed to a greater extent 
to the yield of ions at the wavelength used in these calculations. While differences do exist 
between the and F2 ion yields, we see that the differences are no more than a factor of three. 
We conclude that single ionization of F2 is not suppressed with respect to the single ionization of 
N2. 

C. Single ionization yields from and ^©2 

One explanation postulated for the nonatomiclike single ionization in O2 compared to single 
ionization of N2 is the fact that the ground state of O2 is a triplet state with a half-filled open- 
shell structure whereas the ground state of N2 is a singlet state with a closed-shell structure lO]. 
Miith-Bohm et al Isl argued that such a difference in the spin states is unlikely to influence the 
ionization signal since the spin-degrees of freedom are not effectively coupled to a dipole field. 
We are unaware of any calculation to date in which single ionization of O2 is studied as a function 
of the multiplicity of the ground state. 

Table Hill compares the single ionization yields of ^02 and ^62 at the end of a 24 cycle laser 
pulse (r = 31. 2fs) having a wavelength of A = 390nm over the intensity range / = 1 x 10'"^ W/cm^ 
to / = 8 X 10^^ WW. We see that no significant difference of the yields due to the two config- 
urations is evident. Hence, we conclude that the suppression of O2 with respect to that ofXe^ is 
not due to the ground state 0/O2 being a triplet state. 

D. Orbital response of N2, O2 and F2 

In order to gain a better understanding of the ionization dynamics we have investigated the 
time evolution of the Kohn-Sham orbitals. It must be stressed that these orbitals do not have any 
physical significance with the molecular orbitals of the actual system. However, studying the evo- 
lution of the Kohn-Sham orbitals allows us to obtain information about the orbital symmetries. In 
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figure [U we present results for the time-evolution of the Kohn-Sham orbitals of N2 during its inter- 
action with a 24 cycle laser pulse of wavelength A = 390nm. Two laser intensities are presented: 
/ = 1 X 10''* W/cm^ and / = 6 X 10'^ w/cm^. It can clearly be seen that the orbital having the 
same symmetry (Sag) as the valence orbital of N2 shows the dominant response to the field. The 
same behaviour is observed for all other laser intensities considered in table |ll| Thus it is apparent 
that single ionization of N2 occurs predominantly by ionization of the valence electron. Since the 
valence orbital of N2 has a bonding character then, based upon the symmetry arguments of IMST, 
we would expect no ionization suppression to occur - as is the case. 

In figure |2l we present results for the time-evolution of the Kohn-Sham orbitals of O2 during 
its interaction with a 24 cycle laser pulse of wavelength A = 390nm. Again, two laser intensities 
are presented: / = 1 x 10^'^ W/cm^ and / = 6 x lO'"* W/cm^. As in the case of N2 we see 
that the orbital having the same symmetry (l^g) as the valence orbital of O2 shows the dominant 
response to the field and thus single ionization of O2 occurs predominantly by ionization of the 
valence electron. Since the valence orbital of O2 has an anti-bonding character then, based upon 
the symmetry argument of IMST, we would expect ionization suppression to occur - again, in 
accordance with the numerical results and experiment. 

In figure |3] we present results for the time-evolution of the Kohn-Sham orbitals of F2 during 
its interaction with a 24 cycle laser pulse of wavelength A = 390nm. Four laser intensities are 
presented: / = 1 x 10'^ W/cm^ / = 2 x 10'^ W/cm^, / = 4 x lO'^ W/cm^ and / = 6 x 10'^ W/cml 
It is apparent that the orbital having the same symmetry (Itt^) as the valence orbital of F2 does not 
in this case show the dominant response to the field. Instead we see that for all of the intensities, 
apart from / = 1 x 10'"* W/cm^, the orbital having the symmetry Sag shows the dominant response. 
Since this orbital has a bonding character then, based upon the predictions of IMST, ionization 
suppression will not occur. We conclude that single ionization of F2 does not occur predominantly 
by ionization of the valence electron. Thus, molecular structure effects for this molecule are of 
crucial importance. 

In understanding why the orbital having the same symmetry (Ivr^) as the valence orbital of F2 
does not respond predominantly to the field we consider the response of F2 at the laser intensity 
of/ = lxlO''* W/cm^ where the Sag orbital did not show the dominant response to the field. In 
this case the Ivr^ anti-bonding orbital and the Itt^ bonding orbital showed the dominant response. 
Indeed we see that over all laser intensities these two orbitals response in a similar fashion to the 
field. It would therefore appear that an interference is occurring between these two orbitals. 
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To gain further insight we have repeated the TDDFT calculations at the laser intensity of / = 
2 X 10''* W/cm^, whereby only a subset of the orbitals were allowed to respond to the field, all 
other orbitals being frozen. The results are presented in figure IH We see from figure |4tc) that 
when only the l7r„ and lUg orbitals respond to the field, ionization is suppressed. In figures Eta) 
andEtb) we see that the ionization of either the Itt^ or lug is greater than that of the Sag orbital. 
In figure Iltd), when the Sag, Ivr^ and Ivr^ orbitals respond, ionization of the Sag is dominant. 
This behavior illustrates a rather complicated correlated response of the orbitals to the laser field. 
Moreover, from figure El we can see that a resonance effect is coming into play. This is evidenced 
at the laser intensities of / = 4 x 10^"^ W/cm^ and 1 = 6x10''* W/cm^ where the rate of depletion 
of the Ittu and lUg orbitals is not exponential but instead shows evidence of two distinct rates 
being present, namely an intermediate resonance state is initially populated and then undergoes 
ionization itself. Figure |5] presents results for the time-evolution of the Kohn-Sham orbitals of F2 
during its interaction with a 24 cycle laser pulse of wavelength A = 300nm at laser intensities 
7=4x10''* W/cm^ and 1=6x10''* W/cm^. In this case the orbital populations are depleted 
exponentially, thus confirming the presence of the resonance at A = 390nm. 

Hence, we confirm the connection of the suppression of ionization with destructive interference 
of outgoing electron waves from the ionized electron orbital as put forward first by Muth-Bohm 
et al. [8]. However, the prediction of ionization suppression justified within the IMST approach 
through the symmetry of the HOMO is not reliable, since it turns out that, e.g., in the case of F2 the 
electronic response to the laser pulse is rather complicated and does not lead to dominant deple- 
tion of the HOMO. Therefore, the symmetry of the HOMO is not sufficient to predict ionization 
suppression. However, at least for F2, the symmetry of the dominantly ionized orbital is consistent 
with the non-suppression of ionization. 

IV. SUMMARY 

Single ionization of N2, O2 and F2 by intense laser fields having a wavelength of A = 390nm 
has been simulated using a time-dependent density functional approach within the exchange-only 
local density approximation. The results obtained for the single ionization of N2 and O2 are in 
agreement with IMST results [8], if the HOMO is dominantly ionized. This is the case for N2 and 
O2, where for the latter molecule destructive interference leads to a suppression of ionization. In 
addition it is found that the suppressed ionization of O2 with respect to its companion noble gas 
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atom (xenon) does not arise from the multiplicity of the ground state of O2. 

For single ionization of F2, it is found that the ionization signal is not suppressed with respect 
to the yield, in keeping with other static-field TDDFT results lisl . Multielectron molecular 
correlation leads to the dominant ionization of an orbital with symmetry different from the one of 

n 

the HOMO. This is the reason why the prediction based on the symmetry of the HOMO fails BSD. 

The present calculations are carried out within xLDA. Implementations of TDDFT which go 
beyond the xLDA approximation to a self-interaction free exchange-correlation potential improve 
the description of the ionization potentials and will be the focus of future work. In addition only 
parallel transition have been considered in this paper which is a reasonable approximation for 
strong field ionization of diatomics. Future work will focus on generalizing the calculations of 
single ionziation rates to arbitrary orientations between the molecular axis and the laser polariza- 
tion direction. 
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Ionization Potential (eV) 



N2 O2 F2 
Present 15.91 11.45 14.14 
Exact 15.58 12.07 15.69 



TABLE I: Ionization potentials of N2, O2 and F2 calculated using a TDDFT approach within the exchange- 
only local density approximation compared with the experimental values I23]. While the ionization poten- 
tials of N2 and O2 agree to within 5%, those for F2 differ by 10%. 



Laser Intensity 


Ion Yields 


( X lO^^^ W/cm^) 


N+ 


F+ 


1.0 


0.2479176 


0.1051454 


2.0 


0.3804471 


0.2523075 


4.0 


0.4025715 


0.3742050 


6.0 


0.3501940 


0.3278126 


8.0 


0.2919808 


0.3375370 



TABLE II: Ion yields of N2 and F2 after interaction of the neutral molecules with a 24 cycle laser pulse 
having a wavelength of A = 390nm for a range of laser intensities. The ion yields differ by factor of less 
than 3. Hence, the F^ yield is not suppressed with respect to the yield. 
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Laser Intensity Ion Yield 



1.0 


0.0009385 


0.0006335 


2.0 


0.0022583 


0.0017794 


4.0 


0.0038756 


0.0032043 


6.0 


0.0041731 


0.0033902 


8.0 


0.0036198 


0.0027161 



TABLE rH: Ion yields of (assuming either an initial singlet or triplet configuration of the neutral O2 
molecule) after interaction with a 24 cycle laser pulse having a wavelength of A = 390nm for a range of 
laser intensities. The ion yields roughly agree, indicating that the multiplicity of the ground state of O2 do 
not play a role in the suppression of the yield. 
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FIG. 1: Kohn-Sham orbital populations for N2 during interaction with a 24 cycle laser pulse having a 
wavelength of A = 390nm and laser intensity (a) / = 1 x lO^"^ W/cm^ and (b) / = 6 x 10^"^ W/cm^. For all 
laser intensities it is seen that the Kohn-Sham orbital having the same symmetry (Scjg) as the valence orbital 
of N2 shows the predominant response to the field. 
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FIG. 2: Kohn-Sham orbital populations for O2 during interaction with a 24 cycle laser pulse having a 
wavelength of A = 390nm and laser intensity (a) / = 1 x 10^'^ W/cm^ and (b) / = 6 x 10^"^ W/cm^. For all 
laser intensities it is seen that the Kohn-Sham orbital having the same symmetry (1%) as the valence orbital 
of O2 shows the predominant response to the field. 
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FIG. 3: Kohn-Sham orbital populations for F2 during interaction with a 24 cycle laser pulse having a 
wavelength of A = 390nm and laser intensity (a) / = 1 x lO^^^ W/cm^, (b) / = 2 x lO^^^ W/cm^, (c) / = 
4 X IQi"^ W/cm^ and (d) / = 6 X 10^"^ W/cm^. For all laser intensities, except / = 1 x lO^"^ W/cm^, it is 
seen that the Sag orbital predominantly responds to the field. At the lowest intensity the Ivr^ anti-bonding 
orbital (having the same symmetry as the valence orbital of F2) is equally dominant with the Itt^ bonding 
orbital. These two orbitals respond almost identically at all laser intensities. 
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FIG. 4: Kohn-Sham orbital populations for F2 during interaction with a 24 cycle laser pulse having a 
wavelength of A = 390nm and an intensity of / = 2 x 10^"^ W/cm^. In these calculations only the orbitals 
shown in each figure respond to the field. All other orbitals are frozen. It can clearly be seen in (c) that 
when only the Itt^ and liVg orbitals respond to the field, ionization is suppressed. In (a) and (b) we see that 
the ionization of either the Itt^ or liTg is greater than that of the Sag orbital. When the Sag, ItTu and l-Kg 
orbitals respond, however, ionization of the Sag orbital is dominant. 
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FIG. 5: Kohn-Sham orbital populations for F2 during interaction with a 24 cycle laser pulse having a 
wavelength of A = 300nm and intensity (a) / = 4 x lO^'^ W/cm^ and (b) / = 6 x lO^"* W/cm^. At this 
laser wavelength we see that the the orbital populations fall off exponentially, unlike the response at A = 
390nm observed in figure |3l We conclude that the response at A = 390nni is due to the population of an 
intermediate resonance state. 
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